Dark polariton-solitons in semiconductor micro cavities 
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We report the existence, symmetry breaking and other instabihties of dark polariton-sohtons in 
semiconductor microcavities operating in the strong couphng regime. These half-hght half-matter 
solitons are potential candidates for applications in all-optical signal processing. Their excitation 
time and required pump powers are a few orders of magnitude less than those of their weakly coupled 
light-only counterparts. 



Polaritons are mixed states of photons and material ex- 
citations and are well-known to exist in many condensed 
matter, atomic and optical systems jll, [sl, 13 , 11] • We are 
dealing below with a semiconductor microcavity, where 
polaritons exist due to mixing of quantum well excitons 
and resonant microcavity photons [3|, y, |5|] . In the strong 
coupling regime photons, emitted as a result of electron 
transitions, excite the medium and are re-emitted in a 
cascaded manner, which gives rise to so-called Rabi oscil- 
lations [1, 3, 4]. This phenomenon results in the two peak 
structure of the microcavity absorption spectrum. The 
measured spectral width of the peaks corresponds to the 
picosecond polariton life time [3] . This is in contrast with 
the more usual weak-coupling regime (typical for opera- 
tion of vertical cavity surface emitting lasers (VCSELs) 
0), where the slow (nanosecond) carrier dynamics does 
not catch up with the fast (picosecond) photon decay. 
Thereby most of the photons leave the cavity as soon as 
they are emitted. In this regime the response to a pulse, 
resonating with a cavity mode, results in a single spectral 
peak. Thus any potential application of microcavity po- 
laritons in optical information processing leads to a 2-3 
orders of magnitude response time reduction relative to 
the VCSEL-like operating regimes. 

One of the topics of the recent research into the weakly 
coupled semiconductor microcavities has been the lo- 
calised structures of light or cavity solitons 0, 0, @, 0, 
[iqI , im , 0], which have demonstrated rich physics and 
have been proposed for information processing applica- 
tions 0, S]- In the weak coupling regime formation of 
polaritons is irrelevant, since the dispersion of linear ex- 
citations is purely photonic. Slowness of the light-only 
cavity solitons is an outstanding problem, which can be 
rectified in the strong-coupling regime, where potentially 
much faster, but not yet reported, light-matter solitons 
are expected. 

In the last few years extensive studies of the polari- 
tons in strongly coupled microcavities have been strongly 
motivated by the smallness of the polariton mass lead- 
ing to observation of the polariton Bose-Einstein con- 
densation at few Kelvin temperatures [13, 14]. Polari- 
tons have also been recently observed even at the room 
temperatures, see, e.g., [iS^, which has further boosted 



their potential for practical applications. Another very 
important feature of polaritons in semiconductor micro- 
cavities is their strong repulsive interaction (two-body 
scattering) resulting in a substantial defocusing nonlin- 
earity [J 0] . Amongst nonlinear effects predicted or 
observed with microcavity polaritons are optical bista- 
bilit y [TgI, [itI, fisl, [T^ and parametric conversion 
BllilS, 111- Observation of these effects with 
polaritons requires pump intensities of ~ lOOW/cm^ or 
below (see, e.g.. Fig. 1 in [17]), which is less than 
the typical pump of lOkW/cm^ required for semiconduc- 
tor microcavities operating in the weak-coupling regime 
SaiiEP (see, e.g.. Fig. 5 in [9]). 

Solitonic effects with polaritons in bulk media have at- 
tracted a significant (mostly theoretical) attention since 
70s till now, see, e.g., Ref. [2, 25]. In the latest wave of 
research on exciton-polaritons in strongly coupled micro- 
cavities the solitonic effects have not been much of a focus 
yet, with an important exception of a recent experimen- 
tal paper [26| . In this work the authors claim observation 
of dark and bright localized structures or cavity solitons 
in a strongly coupled semiconductor microcavity. Some 
other papers have reported localisation of microcavity 
polaritons due to linear defects [14, .27] , as a result of 
switching between two polarizations [2q], or neglecting 
such important requisites of passive cavities as losses, ex- 
ternal pump and hence bistability [29]. For studies of spa- 
tially dependent polariton dynamics, see, e.g., [sO]- Our 
work is aimed at filling an existing gap in the theoretical 
knowledge about microcavity polariton-solitons. This is 
necessary not only for backing so far limited experimen- 
tal observations [26] , but also and mainly for guiding the 
future work in this direction. 

The widely accepted dimensionless mean-field model 
for excitons strongly coupled to the circularly polarized 
cavity photons is 



dtE 
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iA + il^'l^)^' = iE. 



(1) 



Here E and ^ are the averages of the photon and ex- 
citon creation or annihilation operators. Normalization 
is such that {VLR/g)\E\^ and {VLR/g)\^\^ are the photon 
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FIG. 1: Polariton dispersion calculated from Eq. (|2]), 
ReA±(kx), and renormalized back into physical units. Pa- 
rameters are n = 3.5, operating wavelength A = 0.85//m, 
7^^o 0.1, HQr 2.5meV. 



and exciton numbers per unit area. Here, Qr is the Rabi 
frequency and g is the exciton-exciton interaction con- 
stant. A = {lo — U0r)/^R describes detuning of the pump 
frequency uj from the identical resonance frequencies of 
excitons and cavity, uOj.- Time t is measured in units of 
I/VLr. 7c and 70 are the cavity and exciton damping 
constants normalized to ^r. Transverse coordinates x, 
y are normalized to the value xq = ^/ c/2knVtR where 
c is the vacuum light velocity, n is the refractive index 
and k = nuj/c is the wave number. The normalized am- 
plitude of the external pump Ep is related to the phys- 
ical incident intensity line as \Ep\^ 



gjJinc/^O^R^ 

[3l| . As a guideline for realistic estimates one can use pa- 
rameters for a microcavity with a single InGaAs/GaAs 
quantum well: HQr 0^ 2.5meV^ hg 0^ lO^^eF/im^, see 
[17, 21, 31]. Assuming the relaxation times of the pho- 
tonic and excitonic fields to be 2.5ps gives 7c, — O-l- 
In accordance with this set of parameters the normal- 
ized driving amplitude \Ep\'^ = 1 physically corresponds 
to the external pump intensity ~ lOkW/cm^ . Optical 
bist ability appears for ~ 0.1, it gives the input inten- 
sity ~ lOOW/cm^. Experimentally the polariton bista- 
bility has been observed for values close or even less than 
lOOW/cm^ EE HE [2^. 

First we briefly summarize important aspects of the 
linear dispersion and bistability properties of the above 
equations. Assuming that ^ ~ ^ik^x+ikyy neglect- 
ing pump and nonlinearity we find the dispersion law of 
cavity polaritons 



k'^ - ^(7c +70) 



(P + i(7o-7c))^ 



(2) 



where k^ = k'^ -\- ky. ReAj^ corresponds to the fre- 
quency of the upper polariton (U-polariton) and ReA- 
to the lower polariton (L-polariton) branch, see Fig. [H 
In the strong coupling regime the gap between U- and 
L-polaritons is greater than the linewidth of the branch 
due to ImA± 7^ 0. 

If E'p 7^ 0, then solitons can exist only on a finite am- 



plitude background (£^(±00) 7^ 0), simply because the 
zero homogeneous solution is absent. Therefore we pro- 
ceed with a brief consideration of spatially homogeneous 
solutions (HSs) and their stability. Then we report the 
existence of various cavity polariton solitons (CPSs) and 
study their stability and instability scenarios. HS having 
bistable dependence from Ep is an important prerequisite 
for the soliton existence. E{Ep) is multivalued provided 
that /(A) > 0, where 



/(A)^A(A2 + 7c'-1) 



y37o(A^+7c' + :^). 



7c 
70' 



(3) 



The cumbersome expressions for the roots of /(A) = 
simplify for 7^ = 70 = and give two bistability inter- 
vals A > 1 and — 1 < A < 0. These two intervals overlap 
with the A intervals allowed by the dispersion relation, 
see Eq. (j2j) and Fig. 1. The bistability in the interval 
— 1 < A < appears because of the nonlinear resonance 
of the pump with the L-polaritons whereas the bista- 
bility in the semi-infinite interval is associated with the 
nonlinear resonance of the pump with the U-polaritons. 
Weakly coupled cavities with defocusing nonlinearities 
exhibit bistability only for A > 0, see, e.g. [l2|. Below we 
focus our attention on the solitons linked to L-polaritons, 
therefore our studies are unique to the strong coupling 
regime. Stability analysis of the HS L-polaritons (A < 0) 
has been previously reported for example in Refs. [2l|. 
The lower state of the L-polariton bistability loop can 
be modulationally unstable within some interval of Ep, 
while the upper state is generally stable, see Figs. [S^a) 
and 3(a). Here, modulational instability (MI) we under- 
stand as the growth of linear perturbations in the form 
^ik^x+ikyV+Kt (^ji^i^ growth rate). As A is chang- 

ing from the bottom of the L-polariton branch towards 
the linear exciton resonance, A = 0, the point of MI is 
moving towards the left edge of the bistability loop and 
finally goes beyond the latter, cf. Figs. 2(a) and 3(a). 

Restricting ourselves to the structures independent on 
the polar angle {6 = arg{x + iy)) we find that the time- 
independent CPSs obey 



d^E ldE\ 
dr'^ r dr J 



(7c - iA)E = £p + (4) 



where r = \J + ^ = ^^/[(7o — ^A) + 
z = is found solving the real cubic equation (7q + 
{z — A)'^)z = \E\'^. z turns out to be a single valued 
function of \E\'^ throughout the range of parameters cor- 
responding to the bistability of L-polaritons. Thus the 
potential problem of ambiguity in choosing a root for z 
is avoided. 

We start our analysis of cavity polariton solitons from 
the case, when the MI point of low state L-polaritons 
is within the bistability interval. In many previously 
studied models bifurcation points of the homogeneous 
solutions have been the sites where localized structures 
branch off [12]. Applying the Newton iterative method 
to Eq. dH) we have found a family of small amplitude 
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FIG. 2: (a) Amplitude of the homogeneous state (HS) 
(black line), max|^(x,^)|^ for bright solitons (blue line) and 
min |^(x, 2/)|^ for dark solitons (red line) shown as functions 
of A = —0.7, 70, c = 0.1. (b) is the zoom of the rectan- 
gular area from (a) showing bifurcations of the dark solitons. 
(c,d) Exciton density distribution |^(x,y = 0)|^ across the 
bright (c) and dark (d) solitons for the points marked by 1, 
2, 3 and 4. Full and dashed lines in (a)-(d) mark stable and 
unstable solutions, respectively. 



bright CPSs emerging from the MI point, see the dashed 
red line in Fig. [2](a). Going towards smaller values of 
Ep^ the CPSs become more intense, see Fig. [2fc). The 
Ep value, at which the lower and upper homogeneous 
states can be connected by a standing ID front, is called 
Maxwell point and this is the point where the branch of 
the bright CPSs terminates {Ep = 0.1748). When the 
pump approaches the Maxwell point the soliton broad- 
ens and its peak intensity tends towards the intensity 
of the upper homogeneous state. We also perform a 
full 2D linear stability analysis of the found structures. 
The linear perturbations are assumed in the general form 
e+(r)e^^^+'"^ + el(r)e-^^^+'"*S where J = 0, 1,2, . . . (32|. 
The resulting Jacobian operator is analysed using finite 
differences in r. The linear stability analysis shows that 
the bright CPSs are unstable with respect to the per- 
turbation with the azimuthal index J = and that the 
development of the instability splits the CPS into 2D 
moving fronts. When Ep is close to the Maxwell point 
this instability is relatively weak and bright CPSs can 
be easily stabilized by the spatial inhomogeneities of the 
pump or cavity detuning. This problem deserves more 
detailed investigation and it will be analyzed elsewhere. 

Because of the defocusing nature of the polaritonic 
nonlinearity dark CPSs, see, e.g., [12], are expected to 
be naturally selected by our system and the instability of 
bright CPSs is not surprising. Dark cavity solitons, have 
been previously studied both theoretically and experi- 
mentally for semiconductor microcavities in the weak- 
coupling regime, see, e.g. 0,[3l,[33]- Unlike fiber solitons, 
the dark cavity solitons have no conceptual disadvantage 
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FIG. 3: a) Amplitude of the homogeneous state (HS) (black 
line) and min |^(x,^)|^ for dark solitons (red and blue lines) 
shown as functions of Ep-. A = —0.5, 70, c = 0.1. Bl and 
B2 mark two branches of dark CPSs. (b) is the zoom of 
the rectangular area from (a) showing bifurcations of the B2 
dark CPSs. (c,d) Exciton density distribution |^(x,2/ — 0)|^ 
across Bl (c) and B2 (d) CPSs for the points marked by 1, 
2, 3 and 4 in panels (a) and (b). Full and dashed lines in 
(a)-(d) mark stable and unstable solutions, respectively. (e,f) 
show development of the symmetry breaking instabilities of 
the CPSs marked as 5 and 6 in (b). 



over the bright ones as information carriers. The branch 
of dark CPSs have been found to detach from the left 
folding point of the bistability loop and tend towards the 
Maxwell point, see Fig.[2](a) and the zoomed area in (b). 
At the onset of their existence the dark solitons are seen 
only as a very deviation from the homogeneous back- 
ground. As Ep tends towards the Maxwell point from 
the right, they become much dipper. Near the Maxwell 
point the dark solitons become very broad and can be 
roughly considered as superpositions of infinitely sepa- 
rated ID fronts (1/r term in Eq. (|4|) can be disregarded 
for large distances and the equation becomes effectively 
ID). It is important to note that the relaxation of the 
fronts towards the upper state happens without oscilla- 
tions, however the relaxation towards the lower state is 
oscillatory, see Fig. [Sfc). Thus pinning of the two fronts 
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and hence stabilization of CPSs is possible only for the 
dark structures (see the thick line in Fig. 3(c)). The 
stable branches of dark CPSs are shown by full lines in 
Fig. [2jb). The unstable ones correspond to the instabil- 
ities with J = 0. 

In the case when the lower branch HS is unstable 
within the whole range of the bistability (A = —0.5) 
we have not found bright solitons, see Fig. 3. This re- 
sult is not surprising because the bright soliton branch 
is expected to detach from the point where the lower HS 
changes its stability. This point is now well out of the 
bistability range, which is another prerequisite for their 
existence. However, we have found two distinct branches 
of dark CPSs marked as Bl and B2 in Figs. [3ja),(b). 
The Bl branch bifurcates subcritically from the folding 
point of the upper homogeneous state. Initially unsta- 
ble (J = 0) CPSs become stable after the turning point. 
Close to this turning point the Bl CPSs have a deep 
like shape, while later they transform into dark rings of 
growing radius Fig. [31(c). Note, that close to the turn- 
ing points additional destabilization of dark CPS hap- 
pens due to linear eigenmodes with complex hi and J = 
(Hopf instability, see Fig.[3Ka)) resulting in the formation 
of oscillating dark CPSs. 

The branch B2 consists of ring shaped structures, see 
Fig. [3^d). The linear stability analysis shows that the 



B2 CPSs can be stable (see the interval marked by 4 
in Fig. [31(b)). However, more often, they are unstable 
with respect to perturbations breaking the radial sym- 
metry, i.e. with J 7^ 0. An example of this instability 
development is shown in Figs.[3l(e,f). The dark ring CPS 
shown in Fig. 3(e) is unstable against linear eigenmode 
with J = 3. The broader CPSs undergo azimuthal insta- 
bilities with larger azimuthal numbers J. For example, 
J = 8 for the concentric ring CPS shown in Fig. [31(f). 

In summary: Following a series of recent experiments 
on observation of microcavity polaritons, we have stud- 
ied the formation of spatially localised polariton-soliton 
structures in the strong coupling regime. In particular, 
our results can be used for the interpretation of the ex- 
perimental measurements reported in [26] , where the po- 
lariton Rabi splitting has been observed simultaneously 
with the formation of various bright and dark localised 
structures. Microcavity polariton solitons reported here 
exhibit a picosecond excitation time and can be observed 
at pump powers few orders of magnitude lower than those 
required in the weak coupling regime of the semiconduc- 
tor microcavities [6", 0, 0, 0, EO, Ell • Thus the light-matter 
polariton solitons have potentially significant advantages 
in all-optical signal processing applications over the light- 
only cavity solitons 0, [O, [l2[ • 
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